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The center-of-mass motion of two two-level atoms coupled to a single damped mode of an electro- 
magnetic resonator is investigated. For the case of one atom being initially excited and the cavity 
mode in the vacuum state it is shown that the atomic time evolution is dominated by the appearance 
of dark states. These states, in which the initial excitation is stored in the internal atomic degrees of 
freedom and the atoms become quantum mechanically entangled, are almost immune against pho- 
ton loss from the cavity. Various properties of the dark states within and beyond the Raman-Nath 
approximation of atom optics are worked out. 

PACS numbers: 42.50.Ct, 42. 50. Fx, 42.50.Vk 



I. INTRODUCTION 

Recent advances in cavity quantum electrodynamics 
have significantly expanded our understanding of the in- 
teraction between matter and the quantized electromag- 
netic field A central topic in these studies is the 
theoretical and experimental investigation of situations 
in which a single atom interacts with a small number 
of modes of the radiation field in high-Q optical or mi- 
crowave resonators. In such a setting, the dynamical be- 
havior of the atom is evidently very different from the 
free-space situation and one can observe phenomena such 
as inhibited and enhanced spontaneous emission (3^] or 
Rabi oscillations between two electromagnetically cou- 
pled states [|] . A natural extension of these studies con- 
cerns the modification of the interaction between two 
atoms in a cavity environment. As the interatomic in- 
teraction is ultimately mediated by the electromagnetic 
field, one can expect drastic effects also in this case. The 
interest in this problem has recently grown, stimulated 
in part by the remarkable experiments of Refs. and 
0. For example, several recent articles have examined 
the mutual coherence of the two atomic dipoles under 
various circumstances [|| [TT| . 

In a further study the modification of the near-resonant 
dipole-dipole interaction between two atoms confined to a 
cavity was investigated in detail fl2|| . As a main result it 
was shown that the familiar concept of the dipole-dipole 
potential ceases to be meaningful under certain circum- 
stances. The purpose of the present paper is to continue 
and extend this work, the emphasis now being put on 
the investigation of the actual dynamical behavior of the 
atoms. In particular, we examine the atomic center-of- 
mass motion under the influence of their interaction with 
the cavity field. In order to work out basic aspects of the 
problem we concentrate here on the model of a short and 
closed optical resonator in which the atoms interact ex- 
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clusively with a single damped standing-wave mode of 
the electromagnetic radiation field. An initially excited 
atom will then spontaneously emit a photon into the cav- 
ity mode and subsequently reabsorb it. Consequently, it 
experiences a random walk in momentum space, i.e. heat- 
ing. Due to photon exchange the atom can also interact 
with and excite its partner in the cavity. These processes 
will cease, of course, as soon as the photon escapes the 
resonator due to cavity losses. 

The analysis of this problem shows that, contrary to 
what one might expect intuitively, the presence of the 
second atom does not simply lead to some quantitative 
modifications in the heating and decay process of the 
first. Rather, it causes qualitative changes in the dynam- 
ical behavior of the system. In particular, one observes a 
tendency of the system to settle into so-called "dark" or 
"quasi-dark" states. These dark states consist of super- 
positions of states in which the initial excitation is stored 
in either atom 1 or atom 2, i.e., entangled states of the 
atoms-cavity system. Due to destructive quantum inter- 
ferences these superpositions are completely — or to a 
large degree — dynamically decoupled from the states in 
which the photon is present in the cavity. Thus they are 
immune — or almost immune — to photon decay. Atoms 
in these dark states can be thought of as a new kind of 
"molecule" largely delocalized and bound by the cavity 
electromagnetic field. The focus of the present article lies 
on an analysis of these dark states, which can be viewed 
as a generalization of the antisymmetric Dickc state of the 
theory of super- and subradiance p3[ . To our knowledge, 
the persistence of the entangled two-atom dark states un- 
der the influence of the atomic center-of-mass motion has 
not been previously discussed in the literature. 

Section II introduces our model and establishes the 
notation. In order to motivate the subsequent analysis, 
Sec. Ill discusses some numerical examples that illustrate 
the role of the dark states and demonstrate their long- 
livedness, even in the case of only approximate darkness. 
Section IV gives a detailed analytical discussion of the 
dark states. We first consider the dynamics of the atomic 
system in the Raman-Nath approximation (RNA) , where 
the atoms are treated as infinitely massive. This allows 
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for a very simple and transparent description of the ef- 
fect. We then remove this approximation and demon- 
strate that certain RNA dark states do remain dark in 
the exact analysis. The decay rates of the other RNA 
dark states are estimated, and the analytical results com- 
pared to numerical calculations. A central result is that 
even though these states are only approximately dark, 
they still have extremely long lifetimes. This should ren- 
der the existence of the quasi-dark states amenable to 
experimental observation, at least in principle. Finally, 
further remarks and conclusions are given in Sec. V. 

II. MODEL 

Our objective consists in studying the center-of-mass 
motion of two atoms confined by a trapping potential 
and interacting with the electromagnetic field inside a 
high-Q cavity. In order to work out most clearly some 
of the basic physical effects observable in this system we 
investigate in the following an idealized model problem. 
Questions of experimental realizability will be discussed 
in Sec. V. 

We consider the one-dimensional motion of two two- 
level atoms of mass M trapped inside an infinite square- 
well potential V(x) with boundaries at x = and x = L. 
The upper and lower internal atomic states |e) and \g) 
are separated in energy by an amount of hajQ. The atoms 
which are treated as distinguishable are also placed inside 
a short and closed electromagnetic cavity that is aligned 
with the atomic trap along the x axis. We assume the 
cavity characteristics to be such that the atomic interac- 
tion with the cavity field can be described as a coupling to 
a single mode. In particular, spontaneous photon emis- 
sion into directions other than the x axis is disregarded. 
On the other hand, the damping of the relevant cavity 
mode due to its coupling to the electromagnetic vacuum 
outside the resonator is taken into account. Based on 
this description, the Hamiltonian of the system is 

H = H a + H c + H r + H ca + H cr (1) 

where H a , H c and H r are the free Hamiltonians of the 
atoms, the cavity mode and the vacuum modes, respec- 
tively. They are given by 

H c = huj c a\a c , H r = ^ Tiuj^a^. (3) 

Here, pj is the center-of-mass momentum and Xj the 
position of the jth atom along the x-axis. The atomic 
pseudo-spin operators Uj are defined by Oj = \g,j)(e,j\. 
The annihilation operators for the cavity mode and the 
vacuum modes are denoted a c and a M , respectively, the 



mode frequencies are to c and w^. The interaction of the 
cavity mode with the atoms and with the vacuum modes 
are described by the terms H ca and H cr . In the dipole 
and the rotating-wave approximation, they read 

2 

H ca = ^ Hg cos(kxj + <j>){aja c + <Jja\), (4) 

3=1 

H cr = n (9*n a l a f + 9^ c al), (5) 

where g = (fkJc/^SoLc) 1 ^ 2 denotes the atom-cavity cou- 
pling constant with L c the cavity length. For a planar 
cavity the mode profile is cosine-shaped with wavevector 
k. The phase angle <p characterizes the relative position- 
ing between cavity mode and atomic trap. The coupling 
constant between the cavity mode and the /ith vacuum 
mode is denoted g^. 

In discussing the atomic time evolution we will mostly 
be concerned with situations in which the center-of-mass 
wave function is spread out over a region of extension 
Aa; large in comparison to the cavity mode wavelength 
2ir /k but small in comparison to the trap length L. For 
small enough times the existence of the trap walls may 
thus be neglected. Furthermore, it is assumed that the 
initial wave function can be ascribed a well-defined mo- 
mentum (j>m,po2) and that the effects of the (small) mo- 
mentum spread around this initial value may be disre- 
garded. From the form of the atom-field coupling 
it follows that a single-atom state with momentum p is 
only coupled to states with momenta p ± fik. In view 
of our initial condition we thus introduce the notation 
|(ii, mi), (12, TO2), n c , {n^}) that denotes a state where 
atom j has internal state ij and momentum q$j + rrijTik 
with integer rrij. Thereby , qoj = mod(poj, hk), i.e., 
< Qoj < fik. The number of photons in the cavity 
and the vacuum mode "/^" are denoted n c and n^, re- 
spectively. 

In case only one excitation is present in the system and 
within the realm of validity of the above approximations, 
the general expression for the system state vector is thus 
given by 

= E^CWWKe, m), (g, n), 0, {0 M }) 

+C 2 , m ,„(£)|(s',m), (e,n),0, {0 M }) 
+C 3 , m , n (t)\(g,m), (g,n), 1,{0 M }} 

+ CW./xWKff. m )' (5, n), 0, {1 M }}}. (6) 

We now proceed to eliminate the reservoir degrees of free- 
dom in the system equations of motion with the help of 
the Born-Markov approximation. This introduces an ex- 
ponential decay rate k/2 = ^g^ and a frequency shift 
A c in the dynamics of the amplitudes C 3m . n - For the 
following, we incorporate this shift into the detuning A 
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between the atomic resonance and the cavity frequency 
and work in the interaction picture with respect to luq. 
The effective Hamiltonian time evolution of the system 
before the photon escapes the cavity is then determined 
by 



2,m,n 



3,m,n 



with 



+Cl ,m—l,n + c 2 

,771,71+1 + C 2 , m ,n-l) 



[(q i + mhkf + (q Q2 + nhk) 2 ]/(2hM) (10) 



+ C3,m-l,n)) 


(7) 


+ C3,m,n-l)) 


(8) 


+" 2 


(9) 



describing the influence of kinetic energy. From Eqs. (|7|)- 
@ one notices a further selection rule. For example, the 
set of coefficients C\^ m ^ n with m,n both even, are only 
coupled among each other and to C2, m ',n', m',n' odd, 
and C^^m",n", m" odd, n" even. Note also that Eqs. (^)- 
@ can be written independently of the phase angle (b. 
In the following we set A = for convenience. 

Another interesting situation arises if one takes the 
existence of the atomic trap boundaries fully into 
account. In this case it is convenient to expand 
the center-of-mass wave functions in terms of the 
eigenfunctions of the atomic Hamiltonian (^), i.e., 
Zsm^qxi/ L)sm(Trrx2/ L)/ L, q,r > 1, which can be 
thought of as specific superpositions of momentum states 
with opposite wave vectors. In general, the coupling term 
H ac introduces transitions from a single-particle eigen- 
state ipq^ — ^2/ L sm(irqxi / L)\g(e)} to an infinite num- 
ber of other states ip e q \ 9 ^ ■ Simple selection rules follow if 
one has k = Nir/L with N a positive integer and = 0. 
Under these conditions one obtains couplings only be- 
tween the single-atom wave functions 



ib 9/e 

V2N- 



ib e,s 



Vq+N 



q+2N 



(11) 



with 1 < q < N. The coupling coefficients are all equal 
besides the one between ipN-q and tp q which is of the 
same magnitude, but of opposite sign. After suitable 
identifications the equations of motion for the probabil- 
ity amplitudes of the two-atom system can thus be cast 
into a form identical to Eqs. (fj])-(H) apart from this sign 
peculiarity. An important special case in the coupling 
scheme of expression ( pT| ) arises if q = N. Under these 
circumstances the sequence terminates at tf> q , the part to 
the left of it being obsolete. This special case is of par- 
ticular importance in the discussion of exact dark states 
beyond the RNA. 



III. NUMERICAL RESULTS 

In order to set the stage for the two-atom problem, 
let us first take a brief look at its one-atom counterpart. 
With the help of the procedure used to derive Eqs. (0)-(||) 
we can obtain a similar set of equations for the one-atom 
system, 



iCi^ m — Ul m Ci^ m + — (C^m+l + C2,m-l) 



(12) 



iC2,m — ( w m + A — iK,/2)C2 : m + ^(^m+l + Cl,m-l) 

(13) 

where the notations used here are defined in parallel to 
those for the two-atom case. In particular, we now have 
w m = (<7o + mhk) 2 1 (2HM) . The excited and ground 
state amplitudes are denoted C\ and C 2 , respectively. 
Equations (|l|) and (|l|) are very similar in structure to 
those used in the discussion of near-resonant scattering 
of two-level atoms from a standing- wave laser field jl4| . 
Physically, they describe the atomic momentum spread 
during the interaction with the cavity mode. If we imag- 
ine the standing wave mode as being composed of two 
counterpropagating running waves we see that during 
an emission-absorption cycle the atomic momentum can 
change by an amount of or 2hk. The change depends 
on whether the photon is emitted into and absorbed from 
the same running wave mode or not. Successive cycles 
thus lead to an atomic momentum spread, i.e. heating. 

This is illustrated in Fig. [j], which shows momentum 
distributions P m {r) — |Ci jm (r)| 2 + |C2, m (r)| 2 derived 
from Eqs. ( |l2| ) and ([l3]) as a function of the discrete mo- 
mentum index m and the dimensionless time r = w rec i, 
with bj rec = hk 2 /(2M) being the recoil frequency. These 
distributions illustrate the effective Hamiltonian time 
evolution of the atom before the photon escapes the cav- 
ity, governed by the nonhcrmitian Hamiltonian 



H eff = H « 



ih-a s c a c 



(14) 



H a and H ca referring now to a single two-level atom. The 
initial conditions for the wave function were chosen as 
Ci, m = <5 m , , C 2 , m = 0, and q Q = 0. Figures 0(a), (b) dis- 
play the case of a lossless cavity {k — 0) and a dimension- 
less atom-cavity coupling constant Q — g/2uj rec — 50. In 
Fig. |l|(a), the influence of the kinetic energy term p 2 /2M 
is neglected (the Raman-Nath approximation) and the 
momentum spread grows linearly in time at a rate pro- 
portional to fir. This should be compared to Fig. [j](b), 
which is for the full model including the kinetic energy 
terms. This illustrates the well-known fact that the RNA 
is only valid for short enough times. Due to the in- 
creasing mismatch between the photon energy and the 
atomic energy increment accompanying a photon absorp- 
tion, the width of the momentum distribution eventually 
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stops growing and begins to oscillate. The effects of cav- 
ity damping are illustrated in Fig. |l](c) and (d), which 
again compare the momentum distributions in the RNA 
and the full model, but for a moderate cavity damping 
rate K 1 — K/ui rec — 50, i.e., k'/CI — 0.4. In this case 
the excited state population is damped on a time scale 
approximately given by 4//i'.[] 

We now turn to the two-atom situation, with the goal 
of determining how the previous results are modified 
when we insert a second atom into the cavity. The dra- 
matic changes brought about under these circumstances 
are illustrated in Figs. ||(a)-(d), which show results of the 
numerical integration of Eqs. They depict the 

momentum distribution of the first atom before the pho- 
ton escape, p2\t) = Z)i=i,2,3-n C*,m,n( T )| 2 1 as a Ac- 
tion of to and r in both the RNA and the full model, 
and in the absence or presence of cavity losses. The ini- 
tial conditions were chosen such that both atoms are at 
rest but atom 1 is in the excited state, atom 2 is in the 
ground state and no photon is present in either the cavity 
or the vacuum modes, i.e., 



Ci,m,n(i — 0) — 5i,l5 m fiSn,0 



The study of the momentum distribution in the pres- 
ence of cavity losses [Figs. ||(c) and (d), again with 
k' = 20] also yields a surprising result. One finds again 
that only a small number of modes are significantly pop- 
ulated. But in addition, and in contrast to the one-atom 
case, after an initial transient evolution the total atomic 
population decays only very slowly, i.e., the photon es- 
cape from the cavity is strongly inhibited by the presence 
of a second atom. In fact, the time evolution of the distri- 
bution still bears a strong similarity to the lossless case. 
Furthermore, the RNA yields a good approximation to 
the full model also in the presence of losses. A further 
increase of the cavity damping rate only leads to minor 
changes in the behavior of the momentum distribution. 

A closer look at the long-time behavior is provided in 
Figs. H There, the total probability P = £ m p£ ] of 
finding the excitation in the cavity (curve 1) is shown 
for the RNA (a) and the full model (b). The parameter 
values are chosen as in Figs. |2](c),(d). After a rapid initial 
transient the probability P reaches a constant value in 
the RNA, whereas it still decays slowly in the full model. 
The curves 2 and 3 show the time evolution of |Ci,o,o| 2 + 
( 15 ) |Ci, ,± 2 | 2 + |Ci,± 2 ,o| 2 + |C 2 ,±i, ± i| 2 (i.e., the central and 



and poi — P02 = 0. The atom-cavity coupling is again 
set to il = 50. As a consequence of the selection rules 
mentioned in Sec. II one has for these initial conditions 



p(i) 



the most highly populated side modes)and of |Ci ; o,o| 
alone, respectively. These curves again demonstrate that 
the spread in momentum is strongly suppressed. 



IV. TWO-ATOM DARK STATES 



for m even and 

p(i) = 



E l^2,m,n| 2 



\Q 



3,m.n I 



for to odd. Figures ||(a) and (b) display the case of the 
lossless cavity. One can recognize two main qualitative 
differences from the corresponding Figs. 0(a) and (b). 
First, the momentum distribution no longer spreads sig- 
nificantly: rather, it remains concentrated in the central 
mode (i.e. m = 0) and a small number of side modes. 
The other modes remain almost unpopulated. Second, 
the comparison between the RNA and the full model re- 
sults shows that the influence of the kinetic energy terms 
now is much smaller than in the one-atom case. Con- 
trary to Figs. |l|(a) and (b), for the time considered they 
only lead to some quantitative modifications but not to 
a qualitative change. This property is of course due to 
the concentration of the momentum distribution around 
to = 0. It also indicates that the RNA is a valuable tool 
in the interpretation of the two-atom behavior. 



1 It should be noted that for large cavity damping k' ^> ft/2 
the decay rate of the excited state population goes to zero. 
This stabilization effect, however, is different in nature from 
the two-atom dark states discussed below. 



The results of Figs. || and || indicate that the atomic 
time evolution is characterized by the appearance of dark 
states which have the initial excitation stored in the 
atoms and which are almost immune to cavity damping. 
In this section a detailed analysis of these dark states is 
given. Before turning to the full problem we first work in 
the RNA, which was shown to provide a useful approxi- 
mate description. 



A. Two-atom dark states in the Raman-Nath 
approximation 

In order to investigate the dark states it is convenient 
to work also in the position-space representation. The 
equations of motion for the position-dependent probabil- 
ity amplitudes Ci(x±, x%, t) read 



id = 
iC 2 = 
iC 3 = 



c) 2 



2 



g cos(kxi)Cs, 



~2M \dxf + dx\ ' ( ! 



' 2 

d 2 



777 «3 + 7772 + A - W2 



2M \dx\ 
-g [cos(foci)Ci + cos(kx2)C2 



dx\ 



(16) 
(17) 

(18) 
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In the first special case discussed in Sec. II (i.e., the 
atomic wave packet well localized inside the trap) these 
equations have to be solved in the domain < x\,x% < 
2n/k and the solution must be of the form 



Ci = exp(zp i£i + ipo2X 2 )Ci 



(19) 



with Ci fulfilling periodic boundary conditions. In the 
second case (trap conditions taken fully into account) one 
has to consider solutions with vanishing Dirichlet bound- 
ary conditions in the domain < xi,x 2 < L. 

In the RNA, i.e., after discarding the spatial deriva- 
tives, Eqs. (|l6|)-(|l8|) decouple spatially and can be solved 
immediately. At a given point (x±,x 2 ) they form a ho- 
mogeneous linear 3 x 3-system of ordinary differential 
equations the eigenvalues of which are given by 



Ai =0, 
A 2 , 3 = -ft/4 - iA/2 



(20) 
(21) 



±y ^ + i^j - g 2 [cos 2 (fcxi) + cos 2 (fca;2)]. 

The existence of the eigenvalue Ai whose real part van- 
ishes independently of the values of x\, 22, and k ensures 
that an excitation initially present in the system has a 
finite probability of remaining in it in the limit t — » 00. 
In particular, if the atomic wave function is given at time 
t = 0by 

\ip(xx,x 2 ,0)) = Ai(a;x,a!2)|e,fl',0, {0^}) 
+A 2 (x 1 ,x 2 )\g,e,0,{0»}) 
+A 3 (x 1 ,x 2 )\g,g,l,{0 fl }) (22) 

then the asymptotic state reached by the "atoms + cavity 
mode" system is characterized by the probability ampli- 
tudes (arranged as a column vector in a self-evident way) 

[cos 2 (fc:ri) +cos 2 (fca;2)n :L (23) 

(A\ cos 2 (kx 2 ) — A 2 cos(fca:i) cos(kx 2 ) \ 
—Aicos(kxi)cos(kx 2 ) + A 2 cos 2 (kx 2 ) 
/ 

Note that this state is not normalized, a result of the 
fact that some of the initial excitation has irreversibly 
escaped from the cavity into the reservoir. 

Expression (23) shows that the asymptotic state does 
not have a contribution from the initial amplitude A3, 
furthermore, the final amplitude in the third channel 
where the photon is present in the cavity vanishes. On 
the other hand, if a state has nonvanishing contributions 
Ai or A 2 it will always evolve into a dark state unless 
A\ cos(kx 2 ) — A 2 cos(kxi). The time scale to reach the 
dark state is determined by the eigenvalues A2 and A3. 

From Eqs. ©-([[I) or Eq. @ it follows that a given 
state is a dark state if and only if it is of the form 



A(xi,x 2 ) 



cos(kx 2 ) 
— cos(kxij 




(24) 



and, in addition, it fulfills the appropriate boundary con- 
ditions. The state (|i|) can be viewed as a generalization 
of the dark state in the Dicke theory of sub- and super- 
radiance [fl3|| . 

In the following discussion we concentrate on the case 
of localized atoms in the sense of Sec. II. If one substitutes 
for the function A of expression j24|) the set of plane 
waves exp(iqoiXi + iqa 2 x 2 ) enp[imkxi + i(n + l)kx 2 ], one 
obtains a family of dark states {|<i m n)} which have a 
simple structure in momentum space, i.e., 

\d mn ) =5(|(e,m),(5,n)) + |(e,m),(ff,n + 2)> (25) 
-\(g,m+ 1), (e, n + 1)) - \(g, m - 1), (e, n + 1)}), 

where we have omitted the occupation numbers of the 
photon modes in the notation of the ket vectors for sim- 
plicity. The dark states \d m n) are truly entangled states. 
Since all permissible functions A can be expanded onto 
the indicated set of plane waves the family {|d mn }} forms 
a basis of the "dark" subspace of the total Hilbert space. 
However, this is not an orthogonal basis as a given \d mn ) 
has a nonvanishing scalar product with four other \d m 'n' ) ■ 
Of particular interest in our context is the question 

of how to characterize the asymptotic state \Dmn 9 ^ e ) 
associated with a given initial state \(e/g,m),(g/e,n)). 
Its coordinate representation can be inferred immediately 
from Eq. ( p3|) , but further insight into the nature of the 
state can be obtained from its momentum distribution. 
Equations show that it is sufficient to study this 

question for the state |-Dqo)> sm ce the distributions for 
the other states can be obtained by a suitable shift of 
indices. In coordinate space the state |-Dqq) is represented 
by 



[cos 2 (fca;i) + cos 2 (kx 2 )] 

x (cos 2 (fca;2), — cos(fca;i) cos(fccc2), 0) T . 

Its momentum-space amplitudes 

Ci/ 2 ,m,n = ((e/g,m), (g/e,n)\DH) 
are determined by 



c l , m , n = —J J a dx 1 dx 2 e~ l(mkxi+nkx ^ 



(26) 



fi(xi,x 2 ) 



cos 2 (kxi) + cos 2 (kx 2 ) 



(27) 



with /1 = cos 2 (fca;2) and f 2 — — cos(kxi) cos(kx 2 ). As 
discussed in Sec. II, c 1( -2), miJl =/= only for to, n both even 
(odd). Evaluating the integrals ( |27| ) one finds that the 
amplitudes ci^ 2m ,o, to > 0, are given by 



Cl,2m,0 — ^m.O + 7T^ m + ^rn-l) 
Z7T 



(28) 



where the numbers I m satisfy the recurrence relation 
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Im = -[(-l) m_1 4i - (6m - 3)/ m _i - (m - l)/ m _ 2 ] 
m 

(29) 

and Iq — 7_i = i7r/2. Further relations between the 
amplitudes Ci,m,n are given by 

+ Cl + C 2 , m +1,„+1 + C2, m +1,„-1 — 0, (30) 
Cl,m,n + Ci jmi „+2 — C2, m +l,n+l — C2,m— l,n+l 

= <Sm,o(<S«,0 + 5fi,-2), (31) 

Ci,m,n — Q,±m,±n (32) 

with m,n both even in Eqs. ( |30| ) and (|3l|). Equation 
( p0|) is a direct consequence of Eq. (^|) whereas Eq. ( pil] ) 
follows from the relation 

Irfoo) - + TO - \D\l) l^ 9 -i,i))/2. 

With the help of Eqs. (p8[)-(p2|) all amplitudes Ci^ m , n can 
be calculated iteratively. In this way, one obtains for 
example 

Cl,0,0 = V 2 ) 

C2,±i,±i = 1/tt - 1/2 « -0.1817, 
ci, ±2 ,o = -Ci,o,±2 = 1/2 - 2/tt « -0.1366. 

An interesting way to determine the scalar products 
(-D^,JZ?qq} with a = eg or ge proceeds as follows [the 
method can also be used to derive Eq. (pi])]. The 
asymptotic state |-Dqq) into which |(e, 0), (<?, 0j) evolves is 
uniquely determined. Any state in the "dark subspace" 
orthogonal to \Dqq) must have vanishing overlap with 
|(e, 0), (g, 0)). If we denote by \Dqq) the state |-Dqo) af- 
ter normalization — remember that the dark state into 
which a given initial state evolves is not normalized — 
we must have that 

\D e °) = \D e °){D e Q °\(e,0),(g,0)). 

Comparing coefficients one obtains that 

(D%\D%) = 0.5, (33) 

i.e., the system has a 50% probability to be trapped in 
that dark state. Using the Gram-Schmidt orthogonaliza- 
tion scheme to construct from |£)^ ) a state orthogonal 
to \Dqq) leads to the conclusion that 

<^nPoO> = <*, m ,n (34) 

with k = 1(2) if a = eg(ge), i.e., the asymptotic dark 
states are non-orthogonal, in general. Equations (|33| ) 
and (Q) can be verified by evaluating the scalar product 
in position space. 

From Eqs. @)-© it can be inferred that 50% of 
the population of the dark state is trapped in the state 
|(e, 0), (#, 0)), while the states \(i,m), (j, n)) with \m\ + 
\n\ < 2 (4) hold 91.3% (96.3%) of the population. This 
observation explains the localization of the momentum 
distributions in Figs. || and |[ 



B. Exact and approximate dark states in the full 
model 

Turning to the full model described by Eqs. (@)-(||) or 
(|l^)-(|l8|), i.e., taking the kinetic energy terms into ac- 
count, it becomes apparent that, in general, the states 
\d rnn ) and \D mn ) are no longer exactly dark. By 'exactly 
dark' we mean being an eigenstate of the full Hamilto- 
nian with a purely real eigenvalue. It is therefore natural 
to ask whether the full model sustains exact dark states 
at all. Interestingly, a complete answer to this question 
can be given for both cases discussed in Sec. II, i.e., for 
atoms localized well inside the trap and for atoms expe- 
riencing the trap boundaries. In the first situation there 
are precisely two exact dark states, which are given by 

\D\) = | do, -i) = (cos(fca;2), — cos(fcxi), 0) T (35) 

and 

I-D2) = sin(fca;i) sin(fca;2)(cos(fca;2), — cos(fcxi), 0) T 
= |d_i,o) - |di,o) + \di,-2) - M-1.-2). (36) 

Dark states thus appear only if the atomic momenta in- 
volved are integer multiples of %k, i.e., if qoi = <?02 = 0. 
For the second case, in which the atomic wave functions 
extend over the whole length of the trap, it can be shown 
that exact dark states can only exist if in the cavity mode 
function of Eq. k = ttN/L with integer N > 1 and 
= 0. Under these conditions there is precisely one such 
state which, in the coordinate representation, is given by 
the first line of Eq. @. 

For a proof of uniqueness of these dark states one can 
start from the observation that also in the full model ex- 
act dark states have to be of the form (|2~4| ) . Additionally, 
they now also must be eigenfunctions of {p\ +p\)/2M 
under the appropriate boundary conditions. One then 
expands both A{x\,X2) and A(x%, X2) cos(fcEi/ 2 ) onto a 
suitable set of eigenfunctions. The fact that in the ex- 
pansion of A(xi,X2) cos(fca; 1 /2) there should only appear 
terms of the same energy imposes severe restrictions on 
the possible forms for the expansion of A(x\,X2). These 
requirements can only be met in the cases indicated. For 
the situation in which the atoms extend over the whole 
trap the breakoff of the coupling scheme ( pd| ) if q — N 
(as outlined at the end of Sec. II) turns out to be crucial 
for the existence of the dark state. 

These considerations imply that most dark states 
found in the RNA become unstable in the full model 
since they are orthogonal to the exact dark states, in 
general. The numerical results of Sec. Ill suggest, how- 
ever, that the corresponding lifetimes are still very long 
so that these states may be regarded as "quasi-dark." 
The examples shown referred to cases in which f2, k' 3> 1 
which is the relevant situation in practice as discussed in 
Sec. V. Under these conditions one may treat the kinetic 
energy term (pf +p|)/2M as a small perturbation to the 
RNA Hamiltonian. Applying standard perturbation the- 
ory one obtains an imaginary correction to the RNA dark 
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state eigenenergies only in second order, which already 
indicates that these states will be long-lived. A crude 
estimate of the second-order imaginary part shows that 
the state \D(d) mn ) acquires a finite decay rate that is of 
the order of 

r m „ ~ uj rec {m 2 + h 2 ) 2 n'/fl 2 . (37) 

Thereby, rh and n have to be understood as typical values 
of m and n appearing in the expansion into center-of- 
mass momentum states. The estimate ( |37| ) assumes that 
k' is not too large in comparison to Q so that the square 
root in expression ( pl| ) is essentially imaginary. 

Hence, consistently with the numerical calculations, we 
find that the lifetime of the 'quasi-dark states' is long 
compared to w~ e * under the condition fi, k' 3> 1. Further- 
more, our estimate implies that the decay rate increase 
rapidly for increasing m, n. This is as can be expected, 
since under these circumstances the dephasing between 
the different momentum eigenstates becomes faster. The 
dependence on k' and Q suggest that the coupling to the 
decay channel becomes more efficient when k' is increased 
and decreased. Figure f| shows the decay of the dark 
states \d mn ) for various values of (m, n), k' , and fl. Their 
evolution qualitatively confirms the dependence (^) of 
T mn on these parameters. Thereby, curve (a) should be 
compared to curves (b), (c), and (d) as in each one of 
these one relevant parameter is changed in comparison 
to (a). 



V. SUMMARY AND CONCLUSIONS 

In this paper we have investigated the dynamics of two 
two-level atoms coupled to a single damped mode of an 
electromagnetic resonator, including the effects of pho- 
ton recoil. We concentrated on the situation where one 
quantum of excitation is initially present in the system. 
A generic feature of the atomic evolution is the appear- 
ance of dark states. These states, in which the excitation 
is stored in the internal atomic degrees of freedom, are 
almost immune to photon decay from the cavity. When 
in a dark state, the two atoms become quantum mechani- 
cally entangled and form a new kind of "molecule" bound 
by the quantum of excitation that they share. The state 
of the compound system can conveniently be described in 
terms of a superposition of different states of well-defined 
center-of-mass momentum. A remarkable characteris- 
tic feature of the dark states is their small momentum 
spread, as compared e.g. to the one-atom situation. This 
property makes their description in the Raman-Nath ap- 
proximation quite accurate. While most dark states be- 
come only "quasi-dark" when this approximation is re- 
moved, their damping rate remains quite long indeed. 

When considering the possible practical realization of 
these states, an interesting question concerns the influ- 
ence of a non-constant atomic trapping potential on the 



time evolution of the dark states. If the trapping poten- 
tials can be arranged to be equal for ground and excited 
states, then one can still obtain dark states in the RNA 
(for the full model it can be anticipated that exact dark 
states will not exist any longer, in general). If, as is nor- 
mally the case, these potentials differ from each other, 
even the RNA will not support dark states. However, as 
Eqs. (^6|)-(p^|) show, in the vicinity of the line x% = X2 
the decay will be significantly decelerated so that a rem- 
nant of the dark-state effect might still be visible under 
such circumstances. 

Let us conclude with a brief discussion of the experi- 
mental feasibility to observe such two-atom dark states. 
Recent cavity QED experiments in the microwave and 
optical domain are described e.g., in Refs. |]5|, |l5| , |l6| . 
They typically involve a low density atomic beam passed 
through the electromagnetic resonator, a situation that 
can be modeled in terms of the localized wave packet 
description of Sec. II. In these experiments the residual 
spontaneous atomic decay rate 7 in the cavity (due to 
coupling to vacuum modes) is approximately one order 
of magnitude smaller than the cavity Rabi frequency g 
and damping rate k, which are both comparable in mag- 
nitude. A single-mode description is thus adequate and 
our system (once prepared in the initial state) would have 
enough time to coherently evolve into a dark state. Fur- 
thermore, the recoil frequency uj rec is also very small in 
comparison to g and k (typically less than a factor of 
10~ 3 ) so the RNA should provide a very accurate de- 
scription. In an experimental realization a main diffi- 
culty would certainly consist in efficiently preparing the 
initial system state. From this point of view, the optical 
regime does not appear as promising as the microwave 
regime: First, due to the short free-space spontaneous 
lifetime of optical transitions the atoms probably could 
not be prepared in the excited state before they enter the 
cavity. Second, if they are both simultaneously excited 
inside the cavity the probability of coupling to the dark 
state is relatively low. 

An experiment involving a microwave cavity might 
proceed as follows. Diatomic molecules in a low-intensity 
beam are dissociated such that the two fragments are of 
nonvanishing opposite spin. The atoms can thus be sep- 
arated in an inhomogeneous magnetic field. One atomic 
beam is subsequently prepared in the Rydberg ground 
state, the other one in the excited state. Using atom op- 
tical elements the two beams are guided such that they 
intersect each other in the microwave cavity (at a small 
angle). As the molecular dissociation creates atom pairs 
it should be possible to arrange the setup so that both 
partners pass the cavity simultaneously with high prob- 
ability. The experimental parameters should be chosen 
such that a single atom always leaves the cavity in the 
ground state. The signature of the formation of a dark 
state would consist in detecting an appreciable fraction 
of atoms leaving the cavity in the excited state. In or- 
der to obtain more information about the nature of the 
dark state one could for example additionally observe the 
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spatial atomic density distribution. 
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FIG. 3. Time evolution of the total probability P(r) 
to find the excitation in the cavity (curve 1), of 
|Ci, ,o| 2 + |Ci,o,± 2 | 2 + |Ci, ±2 ,oi 2 + |C 2 ,±i,±i| 2 (curve 2), i.e., 
the central mode and the most highly populated side modes, 
and of |Ci,o,o| 2 alone (3). The parameter values are A = 0, 
fi = 50, k = 20. (a) RNA, (b) full model. 



FIG. 4. Total survival probability P(t) for initial states 
\dmn) in the full model under various conditions. Parameter 
values (a) m = n = 0, k' = 20, Q — 100; (b) m = n — 0, 
k' = 100, n = 50; (c) m = n = 0, k' = 20, Q = 25; (d) m = 0, 
n = 2, «' = 20, n = 50. 
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FIG. 1. Time evolution of the single-atom momentum dis- 
tribution P m for parameter values A = 0, Q = 50. Initially, 
the atom is at rest in the excited state, no photons are in the 
cavity and the vacuum, (a) RNA and k' = 0, (b) full model 
and k = 0, (c) RNA but k' = 20, (d) full model and k = 20. 
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